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Integral Root A" 1 *'* (DOI or SADOl) 
Closed with one program parameter 
2k 

0-5 



+ 2 



-39 



If the "Integral Root" subroutine starts at location t, 
then it is entered (with A in the accumulator) by the 
following § 



s 


50 pF 
50 sF 


s + 1 


26 tF 



Negligible for the special cases A = |A | > p x 2 
A table of typical times (in milliseconds) follows 



-59 



P/A 


el 


•2 


♦3 


•5 


.0 


2 


35 


30 


35 


25 


30 


3 


55 


*5 


ho 


ko 


30 


if 


55 


60 


50 


50 


hO 


10 


120 


105 


120 


105 


90 



For large p and small A the times are considerably greater • 

This routine computes the pth root (p, a positive integer, 
2 < p < 1023) of a 39 binary digit real argument A, 
-1 < A < 1. If |A| > 1 - p x 2~ 39 , then + ' (l-2~ 39 ) is 
taken as A"^ P o Another special case is A = 0* in which 
case A ' p =0. Otherwise A 1 ' 5 is found by Newton's 
iteration method, in which 

x Q = 1.2~ 39 



n+1 



= * + Vp LA/x p ~ ) — x ] 

n ' * ■ * n ' n 



- 2 - 

Convergence of x to A ' p is assumed when 
n 

1/P t(Vx n P ' 1 )=x n ] >0. 

If p is even, of course, A must be non-negative and in 
this case the non-negative real pth root is found . At 
the end of the routine the accumulator contains the 
signed pth root of A. 
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ORDER 




I0TES PAGE 1 





JDOK (E2) 
40 F 

S5 IF 




N(0) = Argument A 


1 


L4 l 








42 19L 




Set link address 


2 


00 19F 
42 21L 




Set N(21L) = p x 2~^ 


3 


L3 F 




■ 




32 19L 




Exit if A = 


4 


LI 22L 
40 2F 


From 17 r 


Set x Q = 1 - 2 Jy 


5 


40 IF 




W(l) « X 




L7 F 
L4 21L 
LO 22L 




Exit if |A| > 1 - (p-i) x 2~ 59 


7 


36 18L 
19 3'7F 






8 


40 3F 
LO 21L 


From 12 




9 


32 12L ., 
50 IF 




i 

N(l) = x P " 1 | 
v ' n 


10 


75 2F 
40 IF 




x 


11 


L5 3F 
L4 23L 






12 


26 8L 




1 




L7 F 


From 9 





13 


50 22L 

66 IF 




i 

| 


14 


S5 F 
LO 2F 




Form and test l/p(A/x -x ) 
'^ v ' n n' 


15 


10 39F 
66 21L 
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1 


16 


S5 F 
36 18L 








■ 


17 


IA- 2F 
22 ^L 




"i " Vi 


• 


"'■ 


18 


L5 F 


From IflS* 


Test sign of A 






19 


36 20L 
L3 2F 




Exist vith MR^ = + 


|A| 1/P 




20 


22 F 
LT 2F 


By V 

From 20*, 3* 
18' 










22 19L 










21 


00 F 
OOF 


By 2* 


p x 2~ 3 ^ 






22 


80 F 
00 IF 




- 1 + 2~ 59 






23 


00 F 
00 IF 




2 -39 







